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Domain wall fermions provide a complimentary alternative to traditional lattice fermion approaches. By in- 
troducing an extra dimension, the amount of chiral symmetry present in the lattice theory can be controlled in 
a linear way. This results in improved chiral properties as well as robust topological zero modes. A brief intro- 
duction on the subject and a discussion of chiral properties and applications, such as zero and finite temperature 
QCD, J\f = 1 Super Yang-Mills, and four-fermion theories, is presented. 



1. Introduction 

When a continuum theory is regularized by 
the lattice the original number of fermion species 
is increased by 2 d , where d is the dimension of 
space-time, with a net chirality of zero. This 
is the well-known fermion doubling problem |Q. 
For vector theories like QCD, this problem has 
been traditionally treated at the expense of exact 
chiral symmetry by using Wilson Q or Kogut- 
Susskind || fermions. Chiral symmetry is broken 
for any non-zero lattice spacing a, but is recov- 
ered together with Lorentz symmetry as a — ► 
resulting in the correct target theory in the con- 
tinuum. However, since any numerical simulation 
is done at non-zero a, one must be able to sim- 
ulate at small enough a so that the effect of the 
breaking does not obscure the underlying physics. 
The problem is that the computational cost of de- 
creasing a is large. For example, in the full theory, 
decreasing a by a factor of 2 requires a factor of 
28-io more computations. 

In recent years an alternative lattice fermion 
method, domain wall fermions (DWF), has been 
developed and used in several applications. The 
method was introduced by D.B. Kaplan 
and was further developed by Neuberger and 
Narayanan p] and by Shamir and Furman 
Domain wall fermions are defined by extending 
space-time to five dimensions. A non-zero five 
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dimensional mass too is present. The size of the 
fifth dimension is L s and free boundary condi- 
tions for the fermions arc implemented. As a 
result the two chiral components of the Dirac 
fermion are separated with one chirality bound 
exponentially on one wall and the other on the 
opposite wall. For any lattice spacing the two 
chiralities mix only by an amount that decreases 
exponentially as L s — > 00. At L s = 00 chiral sym- 
metry is exact even at non-zero a. In this way 
the continuum a — > limit has been separated 
from the chiral L s — > 00 limit. Furthermore, the 
computing cost increases only linearly with L s . 
Therefore, this method provides a new practical 
"knob" on the amount of chiral symmetry break- 
ing due to the lattice discretization. 

In order to have additional linear control 
over the "effective" fermion mass m e ff, a four- 
dimensional mass rrif is introduced by explicitly 
mixing the plus chirality from one wall with the 
minus from the other. The parameters L s , too 
and TO/ control the effective fermion mass m e ff. 
In free one flavor theory one has § : 

m cS = m (2 - TO )[to/ + (1 - m ) Ls ], 0(m (2 (1) 

The (1 — too) Ls is a mass induced from the non 
zero overlap of the two chiral components. This 
residual mass becomes zero as L s -> 00. If 
2 < too < 4 on has four flavors and if 4 < too < 6 
six flavors. There is a symmetry around too = 5. 
In the presence of increasing interactions these 
ranges get renormalized and shrink in size. Fur- 
thermore, the decay rate is not just the simple 
— ln(l — toq) of eq. |l| but also depends on the 
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gauge coupling in a complicated way. 

The five dimensional theory contains L s heavy 
species. For small L s these would not affect the 
dynamics. However, as L s is increased they may 
introduce bulk effects. These effects must be 
subtracted out by dividing the five-dimensional 
fermion determinant by the determinant of the 
same operator without domain walls, i.e. with 
anti-periodic boundary condition along the fifth 
direction ||. This can be accomplished by intro- 
ducing Pauli-Villars type fields as in || or with 
the modification Four-dimensional fermion 
fields q, q used in Green's functions are con- 
structed from the s = plus chirality and s = 
L s — 1 minus chirality components of the five 
dimensional fields ||. Gauge fields are intro- 
duced by treating the s direction as an internal 
flavor space m. They are only defined in four- 
dimensions, and they do not change along s. This 
construction is used by today's numerical simula- 
tions and is referred to as domain wall fcrmions. 

Because the gauge fields are constant across the 
fifth direction s, the corresponding transfer ma- 
trix T is independent of s. Therefore, the prod- 
uct of the transfer matrices along s is simply T Ls . 
For L s — > oo this becomes a projection operator 
to the ground state defined by the correspond- 
ing four-dimensional Hamiltonian TL. It turns 
out that TL is related to the hermitian Wilson 
Dirac operator j^Dwi'm- = —mo) where Dyy is 
the Wilson operator. Then the fermion determi- 
nant det D of the five dimensional theory contains 
a factor of the form (&|0). The boundary state \b) 
is independent of the gauge field and has fixed 
filling level while |0) is the ground state of TL and 
its filling level depends on the gauge field. This 
construction is the overlap formalism Q and has 
a powerful consequence. If the filling level of |0) 
is different from that of \b) their overlap is ex- 
actly zero. The fermion determinant det D can 
have exact and robust zero modes |]! The cor- 
responding index is simply the difference of the 
filling levels and is an integer. It turns out that: 



index = (K + - K_)/2 



(2) 



where K± is the number of positive/negative 
eigenvalues of TL. For a finite lattice, TL is a finite 
discrete matrix amenable to numerical analysis. 



Even though on the lattice there is no topol- 
ogy in the strict sense it has been found that for 
smooth gauge field configurations the index the- 
orem is obeyed exactly || and that for configu- 
rations generated by numerical simulations of a 
quantum theory the index theorem is obeyed in a 
statistical sense pTo| , |xT[] . Also, it has been found 
that the zero mode properties are maintained to 
a good degree for relatively small L s | Q . 

Unfortunately, the same property that makes 
DWF so attractive is also the reason for a weak- 
ness. When a configuration changes from one sec- 
tor to another the number of negative eigenvalues 
of TL changes and therefore an eigenvalue must 
cross zero. As a result T develops a unit eigen- 
value. In that case there is no decay along s, the 
two walls do not decouple even as L s — > oo and 
the chiral symmetry can not be restored. The set 
of configurations for which TL has a zero eigen- 
value is of measure zero and therefore this is not 
a problem |]||. However, configurations in their 
vicinity have very slow decoupling and as a result 
large values of L s may be needed making the nu- 
merical simulations expensive. This is a problem 
for large lattice spacings ( for a -1 w 650 MeV 
an L s « 100 may be needed) but it diminishes 
rapidly as a becomes smaller (for a -1 « 2 GeV 
an L s k, 20 is adequate for most problems). 

The subject of these proceedings is to report 
on work done during the past year relating to the 
chiral properties of DWF as well as to their ap- 
plications in QCD and other vector-like theories; 
for reviews from previous years please see ]l3| 

2. How many ways? 

Before continuing it is important to put DWF 
into perspective. During the last few years a 
wealth of new lattice fermions with improved chi- 
ral properties have also been found. One could 
now wonder; How many different methods are 
there that solve the doubling problem for vector 
theories? Here is a list: 

1) Wilson fermions, 1975 @. 

2) Kogut-Susskind fermions, 1975 || 

3) Domain wall fermions, 1992-1994 J 

4) Infinitely many fields, 1992, 1998 Jl4 

5) Overlap fermions, 1992-1994 § 
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6) Neuberger fermions, 1997 |iq| . 

7) Perfect action fermions, 1997 0. 

8) Ginsparg- Wilson fermions, 1982, 1997 @,|||. 

9) Molecule chains, 1999 §o|. 

10) Topological QFT in 5D, 2000 ||. 

These methods are substantially different in ap- 
pearance but they are of course related since they 
all describe the same continuum limit. Further- 
more, even at finite lattice spacing many of these 
methods are intimately connected, some with ex- 
act analytical relations. Methods 3-10 share a 
common characteristic: in some limit of their pa- 
rameters an "infinity" is present allowing exact 
chiral symmetry even at non-zero a. 

One may think that such a proliferation of 
methods is unnecessary. However, further dis- 
coveries or improvements may come by exploit- 
ing these different ways of solving the same prob- 
lem. The following is a quote from R. Feynman's 
"The Character of Physical Law" @: "There- 
fore psychologically we must keep all the theories 
in our heads, and every theoretical physicist who 
is any good knows six or seven different theoret- 
ical representations for exactly the same physics. 
He knows that they are all equivalent, and that 
nobody is ever going to be able to decide which 
one is right at that level, but he keeps them in 
his head, hoping that they will give him different 
ideas for guessing." 

A list of references for works presented in this 
conference on new lattice fermions other than 
DWF is given in Q-Q. The reader should be 
aware that a large body of work on new lattice 
fermion methods is available. 

3. Chiral properties 

The chiral properties of DWF are determined 
by how well the chiral modes are localized on the 
walls. How fast the modes decay away from the 
wall is a dynamical property that depends on the 
other parameters of the theory such as couplings, 
masses, domain wall height, and volume, as well 
as the type of pure gauge action that is used. The 
broad properties of this dependence have been 
known for some time. However, work that was 
done during the past year has brought these prop- 
erties into much sharper focus |pq]-||73|. 



Before proceeding, it is interesting to see how 
the decay rate depends on the coupling constant, 
or equivalently the lattice spacing, in a simple 
case: the two flavor Schwinger model. In that 
model the chiral condensate (qq) must be exactly 
zero if chiral symmetry is exact. A non zero value 
is a direct measure of breaking induced by the 
fermion regulator. In figure 1 (qq) in units of the 
photon mass m 7 is plotted versus L s . The data 
are from a dynamical numerical simulation 
From top to bottom they correspond to decreas- 
ing lattice spacings of ~ 1/6, ~ 1/8, ~ 1/10, and 
~ 1/12. The physical volume is kept fixed. There 
are a few observations that can be made: 
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Figure 1. (qq) in units of m 7 vs. L s from a 
numerical simulation of the dynamical two flavor 
Schwinger model with too = 0.9, rrif = 0.0 |J. 

a) The data is consistent with exponential 
restoration of chiral symmetry. 

b) There is a fast decay rate up-to L s » 10, fol- 
lowed by a slower one for larger L s . The anal- 
ysis in H suggested that the fast decay rate 
is controlled by the fluctuations of the gauge 
field within a given topological sector, while the 
slower decay is associated with topology chang- 
ing configurations. However, this also reflects the 
author's bias and does not exclude a variety of 
other possible decay functions. 

c) Both decay rates become faster as the lattice 
spacing is decreased. This is important since oth- 
erwise numerical simulations with DWF would 
not be very attractive. 

Probing the effects of finite L s in QCD is much 
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more challenging because of spontaneous symme- 
try breaking. Much of the progress during the last 
year is due to inventive new probes that measure 
the residual mixing between the two walls. In a 
low energy effective Lagrangian sense one expects 
that this mixing will be reflected in the presence 
of an additive residual mass term m res such that: 



TOoff = rrif + m YCS (S G , /?, L s ,m ,m f , V) 



(3) 



where Sq is the form of the pure gauge action, 
(3 = Q/g 2 with g the coupling constant and V 
is the four-dimensional volume. The various ap- 
proaches are: 
i) From the pion mass: 



m f + r< cs 



s) = (4) 



ii) From PCAC @ and fourth ref. in jjj 
r _ 1)^,(0,0)) 



(5) 



rcs <£*^0M» 1)^(0,0)) 

where is the pseudo-scalar density, and J^ q is 

constructed from fields at L s /2. 

iii) Using the Gell-Mann-Oaks-Renner (GMOR) 

relation ||,|6^|37j one can fit (qq) vs. m/ to the 

form: 



(to/ + TO ros )x7r = (qq) - fro 



(6) 



where x-a is the pseudo-scalar susceptibility and 
bo is a constant. 

iv) By fitting the eigenvalues A of the 5D hermi- 
tian DWF operator vs. m/ to the form: 

A 2 = ti 2 [A 2 + (m f + 5m) 2 ], to a cs = (5m) (7) 

where n, A and 5m are fit parameters 

v) By measuring the lowest eigenvalue A m i n of the 
5D hcrmitian DWF operator only for configura- 
tions where the index is non-zero 

mL = (A mi „) (8) 

The index is measured using the overlap defini- 
tion and the eigenvalue flow method ||. For a 
non-zero index A m i„ should be exactly zero in the 
L s — > oo limit; for finite L s the deviation from 
zero is a measure of chiral symmetry braking. 

Using these probes the dependence of m rcs on 
the parameters of the theory were studied in |p5| - 
|n| . In the following, a collection of sample re- 
sults from these references is presented. The fig- 
ures are originals from the corresponding papers. 



When the symbols used in the figures are different 
from the ones in this article the correspondence 
will be noted in the figure caption. 



+ DWF Ls=8 
O DWFLs=12 
□ DWFLs=16 



Figure 2. Quenched QCD; m* es vs. to from an 
8 , j3 — 6.0 configuration. From top to bottom 
L s = 8, L s = 12 and L s = 16 @. 

The dependence of TO rcs on to/ is weak. For 
example, this can be clearly seen in (65). The 
dependence on Too is shown, for example, in fig. 
2 ^54j. The value of mo is rather large except 
in a plateau region between 1.5 and 2.0. Also, 
in that region its value decreases faster with in- 
creasing L s . The L s dependence can be seen in 
fig. 3. The effect of chiral symmetry breaking 
decreases with increasing L s . Also, an Iwasaki 
improved gauge action improves chiral symmetry 
substantially. The value of to w at L s = 64 with 
the standard plaquette action is reproduced with 
the Iwasaki action with L s as small as 16. This 
will be further discussed in section ||. 

Perhaps the most important issue about DWF 
has to do with the way the chiral limit is ap- 
proached. The overlap guarantees that the chiral 
limit will be achieved as L s — > oo. However, it is 
possible that the allowed range of too shrinks to 
zero size at strong coupling. In a strong cou- 
pling calculation indicated that DWF lose their 
light state at very strong coupling. Also, in ]71j| 
it was found that in the large too and strong cou- 
pling region DWF lose the light state. At weaker 
coupling the allowed range of too is expected to 
be non-zero. But even then the precise way chiral 
symmetry is restored is not fully understood. Ob- 
viously this is very important because one would 
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like to be able to fit the data using the full func- 
tional form. 



• Wilson m =1.65 8 x32 

♦ Wilson m =1.65 16 3 x32 
> Iwasaki m =1 .40 8 3 x32 
□ Iwasaki m =1.65 8 3 x32 
A Iwasaki m =1 .90 8 3 x32 
< Iwasaki m =1.65 1 6 3 x32 

— 140 MeV 



10 20 30 40 50 60 70 80 

Figure 3. Quenched QCD; m^.(m/ — > 0) vs. L s 
with the plaquette action at (3 — 5.7 (filled) and 
the Iwasaki action (open) at (3 = 2.2827. [ pl|j72]| . 



The work of []11|,|23| indicates that the the spec- 
trum of the 4-dimensional Hermitian Wilson op- 
erator has a non zero density around zero. Al- 
though that density decreases fast with decreas- 
ing lattice spacing it will induce some slow decay 
rates. This was also found in [j7(| where decay 
rates as small as 10~ 3 where observed. The work 
of |^,[55| provides further insights: Results us- 
ing the Wilson plaquette action for a" 1 « 1 GeV 
and a -1 w 2 GeV are shown in fig. 4 while for the 
Iwasaki action and for the same lattice spacings 
in fig. 5. Also, data points from different volumes 
are shown indicating that the volume dependence 
is weak. Fits with exponential decay to zero were 
not good for a" 1 « 1 GeV Wilson, a" 1 » 2 GeV 
Wilson and a -1 « 1 GeV Iwasaki while fits with 
exponential decay to a constant were reasonable. 
For a" 1 w 2 GeV Iwasaki an exponential fit to 
zero was reasonable. 

More data points may be needed in order to 
precisely establish the form of chiral symmetry 
restoration as a function of L s for different lattice 
spacings. It is possible that a more complicated 
function of L s is required |9|]7l|. For example, in 
|72[ several fitting functions were used including 
double exponentials. 
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Figure 4. Quenched QCD; m^ s vs. L s with the 
plaquette action at (3 — 5.65 (a -1 s» 1 GeV) and 
f3 = 6.0 (a" 1 ss 2 GeV) HIH- 

RG improved 
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mj.g S vs. L s with the 
a' 1 w 1 GeV) and 



Figure 5. Quenched QCD 
Iwasaki action at (3 — 2.2 
(3 = 2.6 (a" 1 « 2 GeV) @J5 

The above QCD results are for the quenched 
theory. The behavior of the dynamical theory has 
been studied in the context of thermodynamics. 
Similar results are evident there as can be seen 
from fig. 6. 

The appeal of DWF is in the fact that any de- 
gree of chiral symmetry can be achieved even at 
finite lattice spacing (provided that the allowed 
mo region is not of zero size). However, in prac- 
tice the required value of L s may be large mak- 
ing the calculations difficult. An alternative ap- 
proach was introduced in [fT2[ ■ At finite L s it is 
possible to approach the chiral limit by tuning 
rrif so that rrif +m res = 0. This approach is very 
similar to Wilson fermions. Even the Aoki phase 
is present for rrif + m ICS < |7^]24||35| . The same 
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chiral symmetry breaking terms are present but 
with much smaller coefficients. 




0.03 - 
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Figure 6. Dynamical QCD; m^ cs vs. L s with 
the plaquette action on 8 3 x 4 lattices with (3 = 
5.2, mf = 0.02, and too = 1.9 (diamonds). The 
cross/square are m^ cs /m* cs for 8 3 x 32, [3 = 5.325 
(JVf = 4 crossover point) |62|]6"7j ] . 
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Figure 7. Quenched QCD; vs. m/ with the 
plaquette action on 16 3 x 32 lattices with (3 = 
5.7, L s = 48, and to = 1.65 [[§||. The star 
indicates the —m^ es point. 

In fig. 7 @,^2| the approach of the pion mass 
to zero for quenched QCD is shown vs. to/. This 
approach may be distorted by: order a 2 effects, fi- 
nite L s effects in the ultraviolet, topological near- 
zero modes (due to quenching), finite volume ef- 
fects, and quenched chiral logs. In the quenched 
theory configurations with near zero modes are 
not suppressed by the fermion determinant. Be- 
cause of the improved chiral properties of DWF 
the propagator has a near-pole as to/ — > caus- 



ing the pseudo-scalar pseudo-scalar (PP) correla- 
tor to diverge. This dramatic effect is shown in 
fig. 8 for (qq) |72f| . However, the coefficient of the 
divergent term is expected to vanish like 1 / \fV 
where V is the volume. This is evident in fig. 8 
|72|,^6| (a similar divergence is present above the 
quenched deconfining transition but it does not 
disappear as V is increased |2£|). Therefore, the 
pion mass extracted from the PP correlator for 
small volumes may be polluted. This can be alle- 
viated by extracting the pion mass from the axial 
axial (AA) correlator instead J72[ |. 

Another pathology of quenching is the presence 
of quenched chiral logs. This can be accommo- 
dated by fitting to a function that contains a log 
as in fig. 7. This fit gives m* es = 0.0073(10) 
which is consistent with the value calculated us- 
ing PCAC mf es = 0.0072(9) [[72). This provides 
a self-consistency check for the understanding of 
chiral symmetry breaking with DWF. 
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Figure 8. Quenched QCD; —(qq) vs. to/ with 
the plaquette action on 8 3 x 32 (circles) and 16 3 x 
32 (squares) lattices with (3 — 5.7, L s = 32, and 
to = 1.65 |H|]. 

However, it must be stressed that to/+to 10S = 
does not eliminate all chiral symmetry break- 
ing effects |8q| . Power divergent operators intro- 
duce 0(m ICS /a 2 ) effects that must be subtracted 
out. These effects can be clearly seen in Q 
where / 2 TO 2 r /48(TO/ + TO res ) does not extrapolate 
to (qq)(mf + m res = 0). 

4. Improvements 

In this section improvements to DWF are dis- 
cussed. The improved methods achieve the same 
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degree of chiral symmetry restoration but require 
less computing. Several methods have been devel- 
oped during the last year |l[fi|^j7(J, 
Also, work on locality properties during the past 
year can be found in |78|-|8l|. 

As seen in section [|the chiral properties of 
DWF are significantly better when an Iwasaki 
pure gauge action is used. This is remarkable; a 
change in the pure gauge action results in an im- 
provement in the fcrmion sector. This is probably 
because the Iwasaki action "over-improves" and 
therefore reduces the frequency of index changes 
due to small instantons shrinking below the lat- 
tice spacing and disappearing. As mentioned in 
section [j], index changes result in slow decay rates. 
This is indicated by the work of where it was 
found that index changes at values of mo where 
simulations are done are mainly due to small ob- 
jects of size around the lattice spacing. 

Using the Iwasaki action at a -1 k 2 GeV and 
L s = 20 reduces m res by a factor of 10 [^9|. At 
a -1 « 1 GeV and L s = 20 the improvement is 
not as dramatic, m rcs is reduced by a factor of 2 
6l]@. At a- 1 w 650 MeV, which is the lattice 
spacing around the N t — 4 thermal transition, 
and L s = 24, m ros reduces only by a factor of 
1.2. Therefore, the Iwasaki action is effective at 
smaller a. QCD thermodynamics with N t = 4 
does not benefit much from it but at Nt = 6 
(a -1 w 1 GeV) one would expect a factor of two 
improvement. Similar results have been obtained 
with other improved gauge actions pi| , p8[ . 

If these interpretations are correct one would 
expect that actions used for "cooling" studies 
may be beneficial. For example, one could 
attempt to reduce index changes by "over- 
improving" the Iwasaki action even more. How- 
ever, a study at a' 1 ~ 650 MeV with an Iwasaki 
coefficient c = — 1 instead of c = —0.331 did not 
show any significant further improvements p6[ . 
Also, a study at a -1 « 650 MeV with an action 
that restricts the plaquette to 1 — Tr[J7 p ]/3 < 0.8 
did not show significant further improvements ei- 
ther J2]]]. Nevertheless, these actions at smaller 
lattice spacings may produce better results. 

Another set of improvements has to do directly 
with the fermion sector: 

1) The projection method p^,F^,Fol consists of 



projecting out a few eigenvectors with the slowest 
decay rate. This method works remarkably well 
as can be seen from fig. 9. 

2) An overlap inspired method introduces a dif- 
ferent 5-dimensional operator jf?]] . Although the 
computational requirement is similar to DWF 
this method is better for analytical work. 

3) A perturbation theory inspired method uses 
beyond nearest 4D neighbor action p^ . It has 
better decay rates at weaker couplings. 

4) An algorithm that takes advantage of the ab- 
sence of interactions along s uses multi-grid meth- 
ods to speed-up computations 74 1. 



■;;iI---F: :■= :■ " 11=1 v'- :-■ --f-/.: 



L B =30: o N=0, DWF 
L B =30: o N=20 proj. DWF 



Figure 9. Quenched QCD; R = mfXn/iqq) 
vs. m/ , where Xn is the chiral susceptibility. R 
should be unity in the absence of chiral symme- 
try breaking. The bottom points are from unim- 
proved DWF while the top are from improved 
DWF by projecting out the 20 slowest decaying 
eigenvectors ]66|,p6|. 



5. Quenched QCD 

Several applications of DWF to quenched QCD 
were done during the past year fs^-JOg). Most 
of them have been presented in the correspond- 
ing reviews but they are also mentioned here for 
completeness. 

Calculations of /^k have been made in fi5|,j72"f 
(see also the review p£| ) . For example, in |72f| it 
was found that Jk,k are close to the experimen- 
tal values. Two different methods were used to 
calculate /„.; one used the PP correlator and the 
value of mfL while the other used the AA cor- 
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relator and did not require a value for m rcs . It 
was found that (U)pp/(U)aa = 1.00(10). This 
provides a strong consistency check for the deter- 
mination of m res . 

A calculation of the nucleon to rho mass ra- 
tio JT2] gave tojv/to p — 1.37(5) in the a — > 
limit and for a box with size ~ 1.6 fermi in 
each direction (the experimental value is 1.22). 
Also, the chiral condensate was measured to be 
{qq) = (245(7) MeV) 3 at L s = 16 and {qq) = 
(256(8) MeV) 3 at L s = 24 (the phenomenologi- 
cal estimate is (229(9) MeV) 3 ). 
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Figure 10. Quenched QCD; Non perturbative 
renormalization factors vs. (ap) 2 on a 16 3 x 32 
lattice with (3 — 6.0, L s — 16, too = 1.8. The top 
points correspond to Zll/Z^- The lower points 
correspond to factors that are zero when chiral 
symmetry is exact [BUp9[. 



preli 


ninary 




_ * 


1ST (exp.) 

N (exp.) 








i x 






i 

X 




1 s 5 

i * 






X 






X 

X 

X 


x B,+ (nucleon) 






□ B,- 


~+ 
I 


1 1 1 


o B 2 " 
i i i 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
m: quark mass 

Figure 11. Quenched QCD; Baryon masses vs. 
to/ on a 16 3 x 32 lattice with (3 — 6.0, L s = 16 
and too = 1.8 US 

The renormalization factors relevant to weak 
matrix elements were calculated using non- 



perturbative methods in [ piU89|1 (see also the re- 
view Q). In fig. 10 one can see that the Zll/Z\ 
factor is not zero as expected. The other factors 
in that figure are required to vanish if chiral sym- 
metry was exact. Indeed, they are consistent with 
zero demonstrating the good chiral properties of 
DWF. Perturbative calculations of these factors 
were done in ]34] , |65] |. Also, a method that extends 
the range over which operators relevant to weak 
matrix elements are defined was developed in |8q] 



(also see the review 29 



The mass splitting between the nucleon iV(939) 
and its parity partner iV*(1535) was calculated 
in |93||37j and is shown in fig. 11. The data are 
consistent with the experimental values. 




Figure 12. Quenched QCD; vs. to/ with 
the Iwasaki action on 16 3 x 24 lattices with (3 — 
2.6, L s = 16, and m = 1.8 HQ. 

The value of Bk was computed non- 
perturbatively in the RI scheme and matched to 
the MS-NDR scheme in p5[ . It was found to 
be B K (2 GeV) = 0.538(8) on a 16 3 x 32 lat- 
tice with L s = 16. This is lower than the value 
B K = 0.562(14) obtained by §3), probably be- 
cause the non-perturbative renormalization con- 
stant of [jS5| is lower than the perturbative one 
used in |83| (see also the review [p8|). 

Operators relevant to ef/e have also been cal- 
culated (8|-|6lj8| (see also the review @). As 
an example the operator is shown in fig. 12 
Q . The magnitude of this operator is small. 

6. Dynamical QCD thermodynamics 

The study of the QCD thermal transition can 
obviously benefit a great deal by using lattice 
fermions with better chiral properties. DWF are 
a natural candidate for such studies and have 
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been used over the past few years to study Nt = 4 
thermodynamics. This effort crystallized during 
the past year @^|, |§-|o|. The main find- 
ings are: 

1) A small but non-zero U{1)a symmetry break- 
ing above and close to the transition is present. 
The good zero-mode properties of DWF lent va- 
lidity to this result. 

2) Two degenerate flavor simulations on 16 3 x 
4 lattices with L s = 24 gave a transition 
temperature of T c = f63(4) MeV and = 
427(11) MeV, where the p mass was used to set 
the scale. Simulations with ordered and disor- 
dered starting configurations agreed after ther- 
malization indicating the absence of a first order 
transition for this value of m n . 

3) Three degenerate flavor simulations on f 6 3 x 4 
lattices with L s — 24 |3(| gave a transition 
temperature of T c = 160(3) MeV and = 
402(7) MeV. Again, simulations with ordered 
and disordered initial configurations indicated the 
absence of a first order transition for this value of 
m w . The time history of these simulations at the 
crossover point is shown in fig. 13. 



■ \ 

1.006 P- A 



Figure 13. Dynamical QCD with 3 degenerate 
flavors; (qq) vs. the trajectory number on a 16 3 x 
4 lattice at the crossover point (3 = 5.225 with 
L s = 32, rrif = 0.02 and too = 1.9. The top data 
is from a disordered initial configuration while the 
bottom is from an ordered one H] . 

Clearly, m n is very large. From fig. 6 it was 
estimated that L s « 100 may be needed to make 
m-n f» 200 MeV. On the other hand, the transi- 
tion with Nt — 6 will be at a -1 1 GeV. From 
quenched studies at o" 1 ~ 1 GeV, see for exam- 
ple fig. 3, one would expect that the use of an 



Iwasaki action would be beneficial and that one 
could achieve w 200 MeV with L s = 24. 

A study of the equation of state using DWF 
could be less demanding. One could integrate 
along the path proposed by G. Fleming |9(|. In- 
tegrate from [3 — to a (3 above the transition 
using the quenched theory which exactly corre- 
sponds to to/ — 1 because of the subtraction of 
the heavy fields. Keeping (3 fixed at that value 
one can then integrate along a line of decreasing 
m j using the dynamical theory. Since this is done 
at realtively small nanL s ~ 24 may be enough 
to make the calculation physically relevant. 



N xN 4 free DWF 

i =1 .0, m,=0.0, ^=1 , 




Figure 14. Free QCD; energy density e on the 
lattice over the continuum Stefan-Boltzman en- 
ergy density esB vs. the number of sites along 
the time direction Nt- The spatial dimensions 
have 4N t sites, to/ = 0.0, and to = 1.0 @. 

In order to understand how good DWF are with 
N — 4 or even iV t = 6 the free fermion energy 
density is shown in fig. 14 j96|. Due to the sub- 
traction of the heavy modes DWF perform better 
than Wilson fermions at small N t . 

7. Super Yang-Mills 

The N = 1 super-symmetric Yang-Mills 
(SYM) theory can be simulated on the lattice 
using traditional techniques, very much like in 
QCD. Simulations using Wilson fermions have 
already been performed, see for example [ 103 1 . 
However, these methods as in QCD, require fine 
tuning in order to recover the target theory in the 
continuum. 

DWF offer an alternative litjjtoj and JTo|- 
[ lTc| ]. At to/ = 0, L s -> oo DWF forbid a 
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gluino mass term and because all other symmetry 
breaking operators are irrelevant SUSY is recov- 
ered in the continuum limit without fine-tuning. 
Furthermore, DWF for mt > have a Pfaffian 
(Pf) with definite sign for any gauge field config- 
uration. This is important since in a numerical 
simulation the Pf is used as a Boltzmann weight 
which must be of definite sign. Wilson fermions 
at finite lattice spacing do not have a Pf of definite 
sign and numerical simulations are done with |Pf | 
instead. Since Pf in the a — > limit has a definite 
sign one expects to recover the target theory but 
one may worry about non-analyticities resulting 
from taking the magnitude. 




10 



30 



40 



Figure 15. Dynamical TV = 1 SU(2) super 
Yang-Mills; (xx) vs - L s with m/ = and too = 
1.9. The diamonds are from an 8 4 , (3 — 2.3 lattice 
while the crosses from 4 4 , = 2.1 p09| , pTo |. 



For TV = 1 SYM the index is equal to 2vN c 
where v is the topological charge and N c the num- 
ber of colors. As a result, instantons break the 
U{\) chiral symmetry down to Z 2 N a and an oper- 
ator with 2N C gluino fields acquires a vacuum ex- 
pectation value. An interesting question is what 
happens to the remaining Z 2 n c symmetry. It is 
expected that Z 2 n c may break spontaneously to 
Z 2 J3l|]. However, on a torus there are constant 
field solutions with fractional topological charge 
1/N C pffil . Although for large volumes these solu- 
tions vanish, for small volumes they may play a 
role and induce a VEV even for zero mass. In par- 
ticular, for mV(xx) "C 1 the v = ±1/TV C sectors 
exclusively contribute to (xx) M - 

Numerical simulations of TV = 1 SU(2) SYM 



were done in [ 109 1 . The results for m/ — are 
shown in fig. 15. As can be seen, (xx)(L s — > 00) 
has a non zero VEV which is maintained even 
for rather small volumes. The time histories of 
the 4 4 volume and to/ = had "spikes" which 
are indicative of zero mode effects. These results 
seem to support a non-zero (xx) due to configu- 
rations with fractional topological charge. Frac- 
tional topological charge configurations have al- 
ready been found in the quenched theory [ 108 1 . 



8. Fermion scalar interactions 

All DWF fields across the extra direction inter- 
act the same way with the gauge field. The inter- 
action of DWF with scalar fields was studied in 
|n| and was found to be different. That interac- 
tion takes place only along the link that connects 
the boundaries of the extra direction. This re- 
veals a richness in the way different spin particles 
couple to DWF. Four- fermion models were stud- 
ied using large N techniques and were supported 
by numerical simulations with N=2. It was found 
that the chiral properties of DWF in these mod- 
els are good across a large range of couplings and 
that a phase with parity-flavor broken symmetry 
can develop for negative my if L s is finite. 

9. Conclusions 

For the first time domain wall fermions sep- 
arate the continuum (a — > 0) from the chiral 
(L s — > 00) limits. Since the computing require- 
ment is only linear in L s they provide practical 
control over chiral symmetry. Furthermore, they 
exhibit robust zero modes which become exact at 
the L s — > 00 limit. 

DWF provide a complimentary alternative to 
traditional fermion methods and can shed light 
to different regions of the parameter space. DWF 
have found a large spectrum of applications such 
as: QCD thermodynamics, quenched QCD, Su- 
per Yang-Mills, four-fermion theories, and the 
Schwinger model. Also, there are proposals for 
improving DWF to achieve the same amount of 
chiral symmetry with less computations. 

Finally, it should be noted that DWF are just 
one of the many new lattice fermion methods. 
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9. 



This wealth of approaches can only lead to further 
new discoveries. 
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